Infrared electron modes in light deformed clusters 
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Infrared quadrupole modes (IRQM) of the valence electrons in light deformed sodium clusters are 
studied by means of the time-dependent local-density approximation (TDLDA). IRQM are classified 
by angular momentum components A/i =20, 21 and 22 whose /x branches are separated by cluster 
deformation. In light clusters with a low spectral density, IRQM are unambiguously related to 
specific electron-hole excitations, thus giving access to the single-electron spectrum near the Fermi 
surface (HOMO-LUMO region). Most of IRQM are determined by cluster deformation and so can 
serve as a sensitive probe of the deformation effects in the mean field. The IRQM branch A/i =21 
is coupled with the magnetic scissors mode, which gives a chance to detect the latter. We discuss 
two-photon processes, Raman scattering (RS), stimulated emission pumping (SEP), and stimulated 
adiabatic Raman passage (STIRAP), as the relevant tools to observe IRQM. A new method to 
detect the IRQM population in clusters is proposed. 



I. INTRODUCTION 

The most prominent electron mode in metal clusters 
is the Mie surface plasmon 0. Being the doorway to 
many structural and dynamical properties, it has been 
much investigated in the past, for extensive summaries 
see e.g. Besides the dominant Mie plasmon, 

many other electron modes, both electric and magnetic, 
have been predicted but not yet observed in clusters, for a 
brief overview see . Similar modes exist in other many- 
fermion systems. In particular, almost all of them were 
identified experimentally in atomic nuclei, for a review 
see 0- 

The analogy with deformed nuclei suggests that de- 
formed clusters should exhibit a family of low-energy 
(infrared) electron modes which are absent in spherical 
systems. Deformed clusters have a partly filled valence 
electron shell and so sustain low-energy one-electron-one- 
hole (left.) excitations inside this shell. Being mainly 
determined by the deformation splitting, the excitation 
energies are rather small and typically lie in the in- 
frared. Since major shells in light clusters cover electron 
levels with the same space parity, the infrared modes 
should have positive parity. Their softness and deforma- 
tion dependence suggest that these should be infrared 
quadrupole modes (IRQM). In axially deformed clusters, 
IRQM are represented by separate A/i =20, 21 and 22 
branches (where A = 2 stands for total angular momen- 
tum and /z its azimuthal component). 

Our present knowledge on infrared electronic excita- 
tions in clusters is very poor. Even basic IRQM prop- 
erties (spectrum, collectivity, evolution with cluster size 
and deformation, responses to external fields) are un- 
known. It is the aim of this paper to investigate main 
features of IRQM from a theoretical perspective, hope- 
fully delivering useful hints for an experimental search. 

IRQM in light and in heavy clusters have a differ- 
ent origin 0. In the present paper, we will consider 



light clusters or, specifically, free deformed singly-charged 
light clusters. Such clusters have essential advantages 
with respect to IRQM. First, they have an extremely 
low spectral density at low excitation energies, a fea- 
ture which facilitates a discrimination of IRQM in two- 
photon processes. The low spectral density also re- 
duces the unwanted level broadening (Landau fragmenta- 
tion, electron-electron collisions). Beams of size-selected 
singly-charged light clusters are readily available. So, 
these clusters are well suited for measurements. Second, 
as is shown below, IRQM in light clusters are dominated 
by one leh component and can be directly and unam- 
biguously interpreted in terms of the electron spectrum 
in the HOMO-LUMO region. So, IRQM can deliver an 
important information on the underlying mean field level 
structure. Since most of IRQM are induced by cluster de- 
formation, they allow to study the evolution of the elec- 
tronic levels with deformation. And third, there is a close 
connection between A/i =21 branch of the IRQM and the 
electronic scissors mode (SM) [1 ITol [Til fl^ . This gives 
a chance to observe the SM, at least indirectly, through 
IRQM. SM is a universal dipole magnetic orbital mode 
peculiar to deformed systems. It was already observed in 
a variety of quantum systems but not yet in clusters. 

Our calculations were performed within a linearized 
time-dependent local density approximation (TDLDA). 
The ionic background is approximated by a soft deformed 
jellium density, for reviews see 0, The method 

goes back to early studies on cluster spectra [l^. It has 
been successfully applied to study the dipole plasmon in 
spherical and deformed clusters. We will con- 

sider here Na clusters as the simplest test case. However, 
IRQM properties worked out below are of a general char- 
acter and so should appear for all other metal clusters 
with pronounced shell structure for valence electrons. 

IRQM are not accessible by one-photon transitions 
which excite exclusively dipole modes. Thus one has to 
use two-photon processes (TPP) where the target state is 
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populated via an intermediate state by two (absorption 
and emission) dipole transitions. The dipole plasmon can 
serve as an intermediate state. 

Experimental data on IRQM are very scarce and lim- 
ited to heavy clusters (e.g. RS measurements of IRQM in 
embedded silver clusters were reported ^^)- Moreover, 
typical two-photon techniques of atomic and molecular 
spectroscopy are rarely used for clusters. This can be 
partly explained by the fact that clusters have some spe- 
cific properties (see discussion in SecEJ, which makes a 
straightforward application of these techniques question- 
able. In particular, the common methods of detection of 
the level population, based on the radiative decay, seem 
not to be appropriate for clusters. In this paper, we 
will examine applicability to clusters of some widespread 
TPP: Raman scattering (RS), stimulated emission pump- 
ing (SEP), or stimulated Raman adiabatic passage (STI- 
RAP). An alternative method for detection of the level 
population in clusters will be proposed. 

Altogether, the present study has two goals: i) to pro- 
vide a first guide on quadrupole electronic excitations in 
light deformed clusters; ii) to explore the access to non- 
dipole electron modes in clusters by modern two-photon 
techniques. 

The paper is outlined as follows. Section ITT1 describes 
the theoretical framework and choice of the most suitable 
clusters. In section lTTTl a general hierarchy of quadrupole 
electron excitations in clusters is presented and the origin 
of IRQM is clarified. In section llVl the IRQM responses 
in one- and two-photon processes are compared and the 
leh nature of IRQM is demonstrated. In section the 
possibility to observe IRQM in two-photon process (RS, 
SEP and STIRAP) is discussed. Conclusions are drawn 
in section lvTl In Appendix A, the dipole matrix elements 
responsible for the coupling between the states involved 
into two-photon process are analyzed. In Appendix B, 
the scissor mode and its connection with the A/i = 21 
branch of the IRQM are briefly discussed. 



II. CALCULATION SCHEME AND CLUSTER 
CHOICE 

A. Basic spectral properties 

The electron cloud is described by the density- 
functional theory at the level of the local-density approx- 
imation (LDA) for the ground state and time-dependent 
LDA (TDLDA) for the excitations, using actually the 
functional of |20|. The ionic background of the cluster 
is approximated by the soft jellium model allowing for 
quadrupole and hexadecapole deformation |2l|. IRQM 
stay in the regime of small amplitudes. We thus em- 
ploy linearized TDLDA, often called the random-phasc- 
approximation (RPA). The actual implementation in ax- 
ial symmetry is explained in [l7| . The reliability of the 
method has been checked in diverse studies of the Mie 
plasmon in spherical 0] and deformed [TtHIsI clusters. 



The choice of the clusters was dictated by the following 
reasons: i) The clusters should be small enough to pos- 
sess a dilute and non-collective IRQM spectrum. Only 
then the spectrum can be resolved and unambiguously 
related to the single-particle levels, ii) Since IRQM are 
mainly induced by cluster deformation, the clusters with 
a strong deformation (both prolate and oblate) are desir- 
able. The simplest case of axial shape is most suitable for 
the analysis, iii) Shape isomers exhibit different single- 
electron spectra 0, 0], which can result in smearing 
out the low-energy spectral lines. The heavier are clus- 
ters, the more isomers they have So, the light 
clusters with a strictly dominant equilibrium shape are 
preferable. Between them, we should choose the clusters 
whose ground state and isomers have the similar (pro- 
late or oblate) shape. Thus we will minimize the blur- 
ring IRQM spectra, iv) The jellium approximation is not 
correct for very small clusters. This establishes a lower 
limit for the cluster size. 

The cluster choice can be preliminary done by review- 
ing the properties of the dipole plasmon which will serve 
as intermediate state in a two-photon process. Besides, 
description of the dipole plasmon allows to estimate ac- 
curacy of the model. Fig. 1 shows RPA results for the 
dipole plasmon in light axially deformed Na clusters: pro- 
late Naj^, Na^ 5 , Na^ 7 and oblate Na^", Na^~ 9 . The shape 
of each cluster is characterized by equilibrium quadrupole 
and hexadecapole deformations, 62 and 64, obtained by 
minimization of the total cluster energy [171 fl8l l2ll| . To 
simulate temperature and other sources of line broaden- 
ing, the photoabsorption cross section is smoothed by a 
Lorentzian |17j . 

As is seen from Fig. 1, the emerging dipole spectra con- 
sist basically of two peaks, the weaker one (A/z =10) cor- 
responds to oscillations along the symmetry axis and the 
stronger one to the orthogonal mode A^t =11 (in fact, two 
identical modes fx = ±1). In most of the clusters, these 
peaks are well separated by the deformation splitting. 
The most heavy sample Na^" 7 shows some Landau frag- 
mentation of the resonance peaks (distribution of the col- 
lective strength between several RPA levels indicated by 
vertical bars). This is the beginning of a trend which per- 
sists towards heavier clusters. The Landau fragmentation 
overrules the deformation splitting already at N ~ 50 
fl8| . Thus one should stay in the region N < 20 where 
the spectra are still dilute and the Landau fragmenta- 
tion is weak. Within that region, one should pick sam- 
ples with a large deformation. The clusters Na^, Na^ 
and Na^" 5 seem to be the best candidates. Besides, their 
isomeric and ground states have the similar shapes p^ . 
oblate in Na^ and prolate in Na^ and Na^ 5 . The anal- 
ysis of IRQM in the next section confirms this cluster 
choice. 

Fig. 1 shows that in general the folded spectra (solid 
lines) agree nicely with the experiment (triangles) which 
signifies the reliability of the method. However, the jel- 
lium plus RPA description, working well for the heavier 
clusters, worsens for lighter ones. In Na 7 and Na^, the 



3 



calculated spectra are redshifted and do not reproduce 
the structure details. The discrepancy is due to the fact 
that light clusters tend to be more compressed due to 
larger surface tension |24f . But we use here a constant 
Wigner-Seitz radius for reasons of simplicity. Besides, 
the jcllium approximation is certainly too rough for the 
description of details in so light clusters. The results for 
the smallest samples may not reach a quantitative level, 
but they are still useful for qualitative consideration as 
will be done here. Fig. 1 shows that even in the lightest 
clusters the jellium TDLDA sufficiently well reproduces 
the average energy, principle gross-structure, and magni- 
tude of the deformation splitting of the resonance. Such 
accuracy suffices for our present survey. 



B. Two photon process 

We consider a two-photon process running via Xfi = 10 
or 11 branches of the dipole plasmon as the intermedi- 
ate states, see Fig. 2. The branches are assumed to be 
well separated by deformation splitting, which is indeed 
the case for strongly deformed light clusters (see Fig. 1). 
Then the two reaction paths can be disentangled by tun- 
ing the photon frequency. This allows to specify and 
monitor the process. For example, if the reaction runs 
only via A/i =10 plasmon branch, then the population 
of the low-energy A/z =22 mode is forbidden and only 
20 and 21 modes remain to be considered. If one of the 
paths is suppressed, (see discussion for Na^ in the next 
section), then the reaction can be tuned to another path. 
The comparison of the path rates delivers important in- 
formation about the structure of the dipole plasmon and 
IRQM. 

The population of the IRQM is approximately calcu- 
lated as a coherent sum of independent two-step pro- 
cesses, each one being a product of dipole photoabsorp- 
tion and photoemission: 

(7 U (2/i 2 i 2 ) = 53<4i T (0-f lAtiiOo^^l/iiii ^2/x 2 i 2 ) 

ii 
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ii 
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Here |l/ziii > and 2/i2i2 > are RPA states of the dipole 
and quadrupole spectra, respectively. The index i\ runs 
over all RPA states with dipole content in the accessible 
energy interval. The photoabsorption (photoemission) 
dipole matrix elements in JJ| define the coupling between 
the ground and dipole (dipole and quadrupole) states in 
two-photon processes. This coupling determines the Rabi 
frequency, the decisive value in TP P. 

Eq. follows from the general expression for RS 
rate given elsewhere (see, e.g. pl|) if one keeps only 
the Stokes term and neglects the interference between 
the neighboring intermediate states. The Stokes term 



alone suffices to illustrate the population of IRQM in two- 
photon processes. Moreover, already the dipole matrix 
elements in (QJ provide a solid ground for the analysis. 
As for the interference, it should be weak in light clus- 
ters. Indeed, in Fig. 1 the dipole plasmon in clusters 
with N < 20 is represented by a few states separated 
by large energy intervals. In any case, the interference 
mainly leads to smoothing of the response, which can be 
taken into account by a reasonable averaging the results. 
In the present paper, we simulate the interference, tem- 
perature and other smearing factors by folding the results 
by a Lorentzian profile with an averaging width A. The 
spectral width is known to increase with the excitation 
energy. So, we use A =0.25 eV for the high-energy dipole 
and quadrupole plasmons (like in our previous RPA cal- 
culations for the photoabsorption spectra 0, 0, Hsj ^ 
and A =0.1 eV for the low-energy IRQM. 

It worth noting that the photoabsorption and photoe- 
mission dipole matrix elements in (1) have essentially dif- 
ferent structure (see Appendix A for more detail). The 
former is determined by the leh part of the dipole op- 
erator er and is generally large. The latter is given by 
lee (electron-electron) and lhh (hole-hole) parts of the 
dipole operator and its value can broadly vary, depending 
on the structure of the dipole and quadrupole states. 



III. IRQM IN THE HIERARCHY OF 
QUADRUPOLE EXCITATIONS 

The spectrum of the three-dimensional harmonic os- 
cillator provides a useful sortingscheme for the valence 
electron levels in metal clusters 13. l2q| . The levels are 
sorted in perfectly degenerate bunches, the major quan- 
tum shells characterized by the principle quantum num- 
ber J\f = 0, 1,2, .... The shells are separated by appre- 
ciable energy gaps and every shell involves only states of 
the same space parity 7r = (—1)^. This oscillator pic- 
ture is well fulfilled in light clusters and provides still a 
good approximation in medium and heavy ones. In ax- 
ially deformed clusters, the levels are characterized by 
Nilsson-Clemenger quantum numbers |26| v = [Afn z A] 
where Af is the major shell as before, n z is the number 
of nodes along the symmetry axis z, and A is the projec- 
tion of the orbital moment onto the same axis. Following 
this scheme, excitations of valence electrons are charac- 
terized by the AJ\f value, the difference in shell for the 
dominant leh jumps. Quadrupole excitations have even 
parity and are collected into the branches A7V = and 
2. The excitations with larger AJ\f are weak and can be 
neglected. 

Figure 3 illustrates the hierarchy of quadrupole exci- 
tations. The lower panels show the unperturbed pho- 
toabsorption calculated for the pure leh states ignor- 
ing the residual interaction. The upper panels give the 
RPA photoabsorption. Let us first look at the strong 
quadrupole resonance appearing at high frequencies in 
the range 2-4 eV. It is blue-shifted by the residual inter- 
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action and the larger the cluster, the stronger the shift 
[2^|. In heavy clusters this resonance is associated with 
the quadrupole plasmon. The resonance is mainly formed 
by E2 transitions over two major shells (AAf = 2). It 
exists in clusters of any shape and exhausts most of 
the quadrupole strength. The resonance is energetically 
very close to the dipole Mie plasmon (in a simple esti- 
mate lo\ — ujpi \J A/(2A + 1) where uj p i is the frequency of 
the volume plasmon £|). In spite of the energy overlap 
with the dominant El plasmon, the E2 plasmon may, in 
principle, be discriminated by means of angular-resolved 
electron energy-loss spectroscopy (AR-EELS) at electron 
scattering angles ~ 6° |scj . 

In the present study, we are interested not in the high- 
energy E2 resonance but in the IRQM which, being di- 
lute and non-collective, can deliver important informa- 
tion about the electron single-particle spectrum near the 
Fermi (=HOMO) level. IRQM are associated with the 
low-energy AAf = branch created by E2 transitions in- 
side the valence shell. Being of the AAf = origin, IRQM 
can exist only in clusters with partly occupied valence 
shell, i.e. in deformed clusters. In Fig. 3, IRQM reside 
at 0.5-1.5 eV. As compared with the E2 resonance, the 
IRQM spectrum is very dilute. It is represented only by 
a few well separated levels. This prevents mixing of left, 
configurations by the residual interaction and creation of 
collective states. The IRQM persist to keep their leh 
nature and so can be easily identified as particular leh 
configurations. As is seen from Fig. 3, the IRQM have 
very weak quadrupole strength in the photoabsorption. 
But they may be accessible in two-photon reactions. 

Fig. 4 shows single-particle levels and leh quadrupole 
transitions inside the valence shells in Na^ , Na^, and 
Na^. Following this scheme, the IRQM excitations in 
Fig. 3 can be unambiguously identified as particular leh 
states. They include {[101] - [110]} 2 i in Na+, {[220] - 
[200]} 20 , {[220] - [211]} 21 , {[220] - [202]} 22 in Na+, and 
{[220] - [200]} 20 , {[211] - [202]} 21 , {[211] - [200]} 21 , 
{[220] - [202] } 22 in Na+ . 

Electric and magnetic modes with the same A w are 
known to mix in deformed systems (see, e.g., pill]'). The 
mixture of electric E21 and magnetic orbital Mil excita- 
tions is especially interesting since itprovides access to 
the orbital Ml scissors mode (SM) [g [l(J EJ The 
properties of this mode are sketched in Appendix B. 



IV. RESULTS AND DISCUSSION 

A variety of IRQM spectral distributions is shown in 
Fig. 5. The first line of the figure contains the quadrupole 
photoabsorption for IRQM states with A/i =20, 21 and 
22. Though IRQM cannot be observed in photoabsorp- 
tion, the latter is useful in any analysis of electron modes 
and so is worth to be considered. The second line ex- 
hibits magnetic dipole strengths for the magnetic dipole 
scissors mode (SM). The next two lines show the two- 
photon populations when the TPP runs separately 



through A/i =10 and 11 branches of the dipole plasmon. 
The intermediate states in the two-photon process in- 
volve all the RPA dipole states from the plasmon region. 

The quadrupole spectra in Fig. 5 can be easily identi- 
fied by leh transitions exhibited in Fig. 4. Most of the 
transitions connect the levels arising due to deformation 
splitting. The corresponding IRQM are determined by 
the deformation and vanish at the spherical shape. 

The first panel of Fig. 5 shows that IRQM are selec- 
tively active in the quadrupole photoabsorption. Some 
modes (A/i =22 in Na^ and A/i =20 and 22 in Na+g) 
have small E2-transition matrix elements and so are sup- 
pressed in this one-photon process. However, as is shown 
below, these modes can be detected in TPP. 

The calculations show that all the IRQM are almost 
pure leh states, i.e. are dominated by one leh compo- 
nent (with the related transition in Fig. 4. The con- 
tributions of the dominant components typically attain 
99 — 100%. Even in A/i =21 states in Na^~ 5 , which most 
deviate from the leh nature, the contribution of the 
leading leh component are 92-93%. This feature is ad- 
ditionally illustrated in Fig. 6 where the unperturbed 
(the dominant leh configuration alone) and RPA IRQM 
strengths are compared. One sees that in both cases 
IRQM spectra and strengths are very similar. Only the 
A/i =21 states in Na^ show considerable redistribution 
of the strength (keeping, nevertheless, predominantly leh 
structure). This case anticipates the involved picture 
for heavier clusters where the collective redistribution of 
strength is much stronger. 

The second panel in Fig. 5 exhibits the photoabsorp- 
tion cross section for the dipole magnetic SM. The one-to- 
one correspondence between E21 and Mil peaks demon- 
strates the intimate connection between these two modes. 
As is discussed in Appendix B, the scissors mode is 
driven by cluster deformation and vanishes in the spher- 
ical case. Two A/i =21 states in Na]^ represent an in- 
structive example. The lower A/z =21 state determined 
by [211] — » [202] transition between the members of the 
deformation multiplct exhibits an appreciable magnetic 
dipole strength and thus should carry a large SM frac- 
tion. Instead, the higher A/i =21 state is determined by 
a [211] — > [200] transition which takes place even in the 
spherical case. The deformation is not crucial here. This 
state favors the E21 field and so can be treated as an 
ordinary quadrupole mode. Altogether, it is seen that 
deformation-induced IRQM provide access to SM. 

The third and fourth panels in Fig. 5 show the TPP 
cross-section for the reaction paths via A/i = 10 and 11 
branches of the dipole plasmon. It is seen that TPP re- 
sponse considerably deviates from the photoabsorption. 
Besides, the population of IRQM crucially depends on 
the reaction path. For example, it is negligible for A/i =21 
states in Na^ (path 10) and Na^ (path 11). Suppression 
of the TPP transfer in these cases is explained by destruc- 
tive effects in the photoemission dipole matrix element 
(see discussion in Appendix A) . This matrix element can 
vary to a large extent, depending on the particular struc- 
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ture of IRQM and dipole intermediate states. So, for the 
efficient population of the particular IRQM, one has to 
choose the most optimal TPP path. The microscopic 
calculations can be used here as a guide. 

Because of the leh nature of IRQM and their strong 
dependence on cluster deformation, they can deliver a 
valuable spectroscopic information on the electron levels 
in the HOMO-LUMO region as well as on the deforma- 
tion effects in the cluster mean field. For example, the 
measurement for A/i =21 state in Nay" can provide the 
energy difference between the Fermi [101] and particle 
[110] levels (see Fig. 4). If to deduce the Fermi level en- 
ergy from the ionization potential data, one finally gets 
the energy of [110] level. Simultaneously we will estimate 
the deformation splitting of the active subshell. The anal- 
ysis will be even more effective, if to combine the TPP 
data with the photoemission and inverse photoemission 
results for the electron spectra. 



mentation in such clusters is very weak, the levels not so 
broad and the lifetimes are much longer. Light clusters 
are strictly dominated by one equilibrium shape and their 
isomeric shapes are close to the ground state one j^]. So, 
we do not expect a strong blurring of the low-energy spec- 
tra. Beams of size-selected singly-charged light clusters 
are readily available. The light deformed singly-charged 
clusters considered in this paper seem to be most suit- 
able for TPP measurements. The cluster temperature ~ 
100 K could be optimal. Then the thermal broadening of 
electron levels is small and, at the same time, the effects 
of the electron-ion coupling are yet smoothed. 

In the next subsections we will consider some typical 
two-photon processes (RS, SEP, STIRAP) and estimate 
their ability to observe IRQM. 



V. EXPERIMENTAL ACCESS TO IRQM 
A. General view 

Two-photon processes are widely used in atomic and 
molecular physics (for a comprehensive discussion see 
[23). Since atomic clusters are similar to molecules, it 
would be natural to use the same reactions in clusters. 
However, applications of TPP to clusters are very scarce. 
This can be partly explained by our poor knowledge on 
non-dipole low-energy electron modes in clusters, partly 
by peculiarities of clusters. In this section we will discuss 
applicability of traditional TPP methods to clusters with 
a particular accent to IRQM. 

As compared with atoms and molecules, atomic clus- 
ters have at least three features essential for TPP. i) 
The dominant decay channel in clusters is usually not 
radiative. Instead, the levels mainly decay through Lan- 
dau fragmentation (dissipation of the collective motion 
through surrounding leh excitation), electron-electron 
collisions, and electron-ion coupling. This property spoils 
the radiation-based methods (typical for atomic and 
molecular measurements) to detect population of levels 
in clusters. In this connection, we will propose an al- 
ternative detection method based on the photoabsorp- 
tion depletion of cluster beam, which is more suitable 
for clusters, ii) Cluster excitations have extremely short 
lifetimes (10-20 fs for the dipole plasmon in medium clus- 
ters) and usually are rather broad (see e.g. Fig. 1 for the 
dipole plasmon). The short lifetimes inhibit the applica- 
tion of adiabatic transfer methods while large widths im- 
pede maintaining the resonance conditions, iii) Clusters 
demonstrate a strong shape isomerism. Cluster beams 
are usually a statistical mix of clusters of a different 
shape. This blurs the measured low-energy electron spec- 
tra. 

Light clusters with their dilute spectrum allow to cir- 
cumvent some of these troubles. Indeed, Landau frag- 



B. Raman scattering 

RS is one of the simplest TPP. In this reaction, a 
dipole laser-induced transition to an intermediate elec- 
tronic level (real or virtual) is followed by the dipole flu- 
orescence to low-energy levels. 

RS measurements of electron infrared modes in clus- 
ters are very rare. We know only one experiment where 
IRQM were observed in heavy silver clusters embedded 
into amorphous silica films (l9| . The dipole plasmon was 
used as a resonant intermediate state. The further mea- 
surements revealed that IRQM are affected by cluster de- 
formation [3(|. These observations show that the radia- 
tive decay of the dipole plasmon is detectable in spite of 
the strong competition with other decay channels. This 
message is very encouraging for application of TPP to 
clusters. 

RS generally assumes that the coupling between the 
intermediate dipole and final quadrupole states is of the 
same order of magnitude as the coupling between the in- 
termediate and ground states. Only then the final state is 
successfully populated. This means that the absorption 
and emission dipole matrix elements defined in i|A3|l and 
(IA5|) should be of the same scale. The calculations show 
that this is indeed the case: in the clusters considered in 
the present paper most of the absorption and emission 
dipole matrix elements lie inside the interval 2-10 eao (in 
atomic units), i.e. are basically of the same order of mag- 
nitude (for exception of some emission matrix elements 
strictly suppressed due to destructive contributions, see 
discussion in Appendix . This means that IRQM in 
light deformed clusters have a chance to be observed in 
the resonant RS (running via the dipole plasmon). 

RS is effective only if the coupling between intermedi- 
ate and desirable target levels is strong enough. As was 
shown for Na^, this is not always the case. To overcome 
this trouble, we should consider the methods with stimu- 
lated emission into a specific target level (SEP, STIRAP). 
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C. Stimulated emission pumping 

Unlike RS, this method exploits two lasers, pump and 
dump The pump laser is responsible for the first 

photoabsorption step. The dump pulse follows the pump 
one with some delay and couples intermediate to target 
states. If the difference between the pump and dump fre- 
quencies is in resonance with the frequency of the target 
state, then the dump radiation stimulates the emission 
to this state. 

SEP enjoys widespread application for atoms and 
molecules but can encounter troubles for clusters where 
SEP methods to detect the population (based on mea- 
surement of spontaneous photoemission or photoelectron 
yield from certain levels excited by probe lasers) can fail 
for short-lived cluster levels with their small radiative 
width. In this connection, we propose a new detection 
method which does not need any probe radiative proce- 
dure. 

The idea is to use depiction of cluster beam caused 
by the photoabsorption. It is known that photoabsorp- 
tion heats the clusters leading to subsequent evaporation 
of atoms. The recoil effect from the evaporated atom 
forces the cluster to leave the beam. Thus the photoab- 
sorption can be measured in terms of depletion. In the 
scheme we propose, the frequency of the pump laser is 
fixed in resonance with the dipole plasmon while the fre- 
quency of the dump laser scans to get the resonance with 
IRQM at up — uje, = loirqm- Then one can detect the 
IRQM population as a dark resonance in the beam de- 
pletion. Indeed if the dump pulse is not in the resonance 
with IRQM, then a usual photoabsorption induced by the 
pump laser and additionally supported by the dump laser 
takes place. In the case of the resonance with IRQM, part 
of the energy will be used not for cluster heating but 
for stimulated radiative decay from the dipole plasmon 
to the given IRQM. So, the photoabsorption (depletion) 
should demonstrate a sharp decrease, kind of a dark res- 
onance. The deeper the dark resonance, the more popu- 
lated is the IRQM. The photoionization channel plays no 
role because the ionization potential in light clusters is 
higher than u)p + u>£>. If the laser intensities are not too 
high, then multi-photon ionization can be also neglected. 
The scheme is simple and suitable for clusters. 

However, even in the best case, SEP can transfer to 
the target level only 25% of the population [3j| . Spon- 
taneous emission from the intermediate level results in 
an additional leaking. So, it is worth to continue our 
discussion and consider more efficient methods. 



D. STIRAP 

Much better results can be obtained with STIRAP [3J 
l35| which promises up to 100% of the population transfer 
from the initial to the target level. So high effectivity is 
provided by the coherent character of the process and 
the principle possibility to avoid the leaking from the 



intermediate levels. Like the previous method, STIRAP 
also implements two lasers: pump laser to couple the 
initial state |0 > and intermediate state |1 > and Stokes 
laser to stimulate the emission from |1 > to the target 
state 2 >. However, STIRAP is more involved because 
it implies the coherent adiabatic transfer. 

STIRAP has three principle requirements: i) two- 
resonance condition tup — lus — ^2 — Wo which allows 
a detuning A = (uji — ujq) — ujp = (u>i —0)2)— uj$ from 
the intermediate state frequency; ii) counterintuitive se- 
quence of the pulses when the Stokes pulse proceeds the 
pump one; 3) adiabatic evolution. 

Under these conditions, one of the time-dependent 
eigenfunctions of the system is a superposition of ini- 
tial and target bare states only: \bo(t) >= co(i)|0 > 
+C2(t)\2 >. It is reduced to |0 > at the beginning of 
the adiabatic evolution and to |2 > at the end. Hence, 
the system finally finds itself in the target state. The 
intermediate state |1 > is not involved to \bo(t) > and so 
is not populated at any time. Thus any leak in the popu- 
lation is avoided and the transfer is complete. The main 
point is to evolve the system adiabatically, keeping it all 
the time in the state \bo(t) >. STIRAP is widely used in 
atomic and molecular spectroscopy. It is rather insensi- 
tive to precise pulse characteristics. Both continuous and 
pulse lasers can be implemented. 

Let's examine the STIRAP requirements for atomic 
clusters, in particular for the IRQM population via the 
dipole plasmon. 

Two-resonance condition. This condition can be obvi- 
ously maintained even for a broad dipole plasmon. How- 
ever, the role of detuning A (which is used as an addi- 
tional tool to prevent the |1 > admixture to \bo(t) >) 
becomes vague. The dipole plasmon vanishes slowly at 
its flanks and so rather large detuning is necessary. This 
can be done even under the restriction A <C ujp ± cos im- 
posed by the rotation wave approximation. At the same 
time, this point was not yet properly investigated and, 
probably, some fraction of |1 > will admix \bo(t) >. Then 
the population will not be complete, although it still can 
remain large. 

Counterintuitive order of pump and Stokes pulses. The 
order when Stokes pulse precedes the pump one is cru- 
cial for the maximal population |34ll35| . Then the Stokes 
pulse prepares the coherent superposition of the interme- 
diate and final states just before arriving the pump pulse 
to stimulate the subsequent desirable transfer. Stokes 
and pump pulses must overlap and the overlapping time 
At determines duration of the adiabatic evolution. This 
time should not be longer than the lifetime of the dipole 
plasmon, i.e. should not exceed 10-100 fs. Though such 
At is extre mely short, it is quite accessible in modern 
experiments [33| . 

Adiabatic passage. This condition is most tough. Fol- 
lowing 0, 133 it can be formulated as 

QAt > 10 (2) 
where 57 = \/£lp + f7| is the average of the pump and 
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Stokes Rabi frequencies. The Rabi frequency is [34| 



n = - 2 - 20 • 108 MM V 7 ^] 5-1 . ( 3 ) 

n V ceo V cm z 

where d is the dipole coupling matrix element in atomic 
units, / is the laser intensity in W/cm 2 , c is the light 
speed, eo is the vacuum dielectric constant. It is easy to 
estimate (for d ~ 5eao and At —10 is) that the condition 
(J2J is maintained for intensities I > 10 12 W/cm 2 . At so 
high intensities, the multiphoton ionization and fragmen- 
tation of clusters should take place 0], which can spoil 
STIRAP. 

However, there is a way to circumvent the trouble. The 
condition was obtained for one intermediate level. 
But the realistic spectrum of the dipole plasmon consists 
of a sequence of dipole levels. In this case, the STIRAP 
condition should be revised. The realistic case is more 
complicated but, at the same time, opens new possibil- 
ities for STIRAP. In particular, one may loose the adi- 
abatic STIRAP condition and thus the requirements for 
the laser intensity. 

A general case of N intermediate states, each with its 
own coupling and detuning was studied in It was 

shown that the trapped adiabatic state \bo(t) > can be 
created only when the ratio between each pump coupling 
and the respective Stokes coupling is the same for all in- 
termediate states. Following our calculations, this condi- 
tion is unrealistic for atomic clusters. However, the softer 
alternative adiabatic requirements can be formulated. In 
particular, in the general case of arbitrary couplings, one 
may tune the pump and Stokes lasers just below all inter- 
mediate states and thus form so called adiabatic-transfcr 
state which also results in a high if not complete popu- 
lation of the target level. Unlike \bo(t) >, the adiabatic- 
transfer state can have admixtures from the intermediate 
states during the evolution period At and so some popu- 
lation leaking is unavoidable. Nevertheless, we have here 
a solid adiabatic transfer with a high population of the 
target state. Finally, one may conclude that STIRAP can 
be used for the efficient population of IRQM in clusters. 



which are widely used in atomic and molecular spec- 
troscopy but not yet for atomic clusters. It was shown 
that, in spite of some peculiarities of clusters (broad reso- 
nances, short level lifetimes, domination of non-radiative 
decay channels), these TPP can be applied to populate 
IRQM. Besides, a new method to detect the population 
of the target cluster states was proposed. TPP mea- 
surements of IRQM can be supplemented by photoemis- 
sion and inverse photoemission experiments delivering 
the similar information. We hope that our analysis will 
stimulate application of the experimental TPP methods 
of atomic and molecular spectroscopy to atomic clusters. 

Non-dipole electron excitations in clusters represent a 
new promising research field. They deliver interesting 
physics and can serve as a robust test for the theory 
which still paid main attention to integral characteristics 
of clusters but not to so fragile patterns as the single- 
electron spectra. It worth noting that the single-electron 
spectra are very sensitive to many cluster features and so 
can be used as an effective tool for investigation of these 
features. 

Though our calculations have been done for sodium 
clusters, the qualitative results we obtained are of a gen- 
eral nature and should be valid for other metal clusters 
as well. In particular, similar deformation-induced IRQM 
are expected for supported clusters where they can serve 
as sensitive indicators of the interface interaction and 
cluster deformation. Besides, IRQM can be important 
for the quantum transport in clusters where they can 
lead to a resonant transmission. 
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APPENDIX A: DIPOLE MATRIX ELEMENTS 



VI. CONCLUSIONS 

We have presented a first exploration of infrared 
quadrupole electron modes (IRQM) in light deformed 
clusters. Most of IRQM are induced by cluster deforma- 
tion and thus can deliver useful information on deforma- 
tion effects, e.g. on the deformation splitting of electron 
levels. Besides, IRQM are about pure one-electron-one- 
hole excitations and so give access to the single-electron 
spectrum. We explained the origin of IRQM and showed 
that they can be easily identified in the dilute spectrum 
of light clusters. 

In the second part, we examined some typical two- 
photon processes (Raman scattering, stimulated emission 
population and stimulated adiabatic Raman passage) 



The dipole and quadrupole electron states are de- 
scribed as RPA modes 

QU = \ E^ef ^ " ^ea{) (Al) 

eft, 

where 

^Pet-^i - f "'\ ■ &-N x , ff (A2) 

are forward and backward amplitudes characterizing con- 
tributions of electron- hole (leh) configurations to the 
mode. Furthermore, Nx^i is the normalization coeffi- 
cient, f J* is the single-particle matrix element of the 
residual two-body interaction, e e /j and cox^i are ener- 
gies of the leh and RPA excitations, respectively. RPA 
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describes equally well both collective and non-collective 
modes. In the latter case (i.e. in the limit \Xfii >—> 

\leh >), only the amplitude \ip e t l \ — > 1 survives while all 
others vanish. 

The photoabsorption matrix element responsible for 
the coupling between the ground and dipole states is de- 
termined by leh part of the dipole operator and reads 

(l»ih\erY lfl \0 >= fl^itir 1 + ^r 1 ) (A3) 

e h 

where f^^ 1 is the matrix element of the dipole transi- 
tion. If to put /g^ 1 w feh^ 1 an d to take into account 
that all the relevant dipole le/i-levels are blue-shifted by 
the dipole interaction, then 

(l//iii|ery lM |0>«2JV lMlil ^ U 2 eh _ > 2 eh (A4) 

e /j e eh IaMi 

and it is easy to see that all leh levels contribute con- 
structively to the transition. 

Instead, the photoemission dipole matrix element con- 
necting the dipole and quadrupole RPA modes is deter- 
mined by lee and lhh parts of the dipole operator and 
has the more complicated structure: 

(2 / u 2 «2|erYi Al |l/iiii) 

ee' h 

+E/w"E^ 242 ^ in +^r^T i )\ ■ (as) 

hh' e 

The terms with backward amplitudes, ~ </>(/>, are omitted 
in l|A5jl since they are usually small. 

The matrix element (|A5(1 becomes simpler when the 
quadrupole state is dominated by a single leh configura- 
tion, say {e^/^}- As is shown in Sec. IIVI this is indeed a 
common case. Then all the quadrupole amplitudes van- 
ish except of IV'gf 2 * 2 ! ~> 1 an d Eq. (|A5|) is reduced to 

({e^llerYiJl^iii} 

= 4E f^lT + E ■ (A6) 

e h 

Now only the dipole amplitudes including the hole (first 
term) or particle (second term) from the pair {^2/12} 
contribute to the transition. The contribution is con- 
siderable only subject to the large dipole amplitude and 
strong matrix element. (The latter take place if the ma- 
trix element fulfills the asymptotic Nilsson selection rules 
[27ll28| : El^t transition between [Mini Z Ai] and [Afjrij Z Aj] 
states is favored if Mi = Mj ± 1, = 7ij Z ± 1 for /.t = 

and rii z — rij z for fi = 1). Under these tough require- 
ments, only a few terms yield large contributions to i|A6() . 
Depending on the structure of the dipole and quadrupole 
states, the contributions can have different signs and thus 
lead to constructive or destructive results. So, the mag- 
nitude of the photoemission matrix element (and thus 
the coupling between the dipole and quadrupole states) 
can be large (like for the photoabsorption) or very small. 



APPENDIX B: SCISSORS MODE 

The SM is a general dynamical phenomenon already 
found or predicted in different finite quantum systems 
(atomic nuclei, metal clusters, quantum dots, dilute 
ultra-cold gases of Bose and Fermi atoms), see the re- 
view All these different systems have two features 
in common: broken spherical symmetry (deformation) 
and a two-component structure. 

Like the quadrupole modes, the SM separates into low- 
energy (AM = 0) and high-energy (AM — 2) branches. 
In this paper we consider only the low-energy branch 
which is mixed with IRQM. Macroscopically, this branch 
is treated as a small-amplitude rotational oscillations of 
a spheroid of valence electrons against a spheroid of the 
ionic background (hence the name scissors mode). Like 
IRQM, the SM is driven by cluster deformation and its 
energy scale is naturally determined by the deformation 
splitting of the electron levels. In axial clusters with 
a quadrupole deformation 82, the energy and magnetic 
strength of the mode are estimated as Ej : 

lo = 2 -^-N- 1 'H 2 eV, B(M1) ~ N^5 2 (Bl) 

' s 

where r s is the Wigner-Seitz radius (in A), N e is the 
number of valence electrons, and /it, is the Bohr mag- 
neton. It is seen that both the energy and strength are 
proportional to the deformation parameter 82 and vanish 
for the spherical shape. 

In axially symmetric systems, the SM is generated by 
the orbital momentum fields, L x and L y , perpendicular 
to the symmetry axis z and, like the quadrupole mode 
\[i = 21, forms the states ^ = 1+ >. SM strongly 
responds to an external magnetic dipole field. Besides, 
it determines the van Vleck paramagnetism in deformed 
clusters 0, El ITTL IT^ . The experimental search of the 
SM in clusters is still in very beginning [lflj . 

There is an intimate connection between the scissors 
and quadrupole E21 modes in deformed clusters. If fact, 
both SM and E21 mode are parts of a general motion 
of multipolarity A/z = 21. To illustrate this point, we 
expand a single-particle electron state v in terms of a 
spherical basis (nLA) 

*„=[Jtf»,A] = E < L RnL(r)Y LA (n) (B2) 
11 L 

and estimate the SM leh matrix element: 

(%\L x \* h ) oc <5 w A p ,a,^ ( B3 ) 
• Y,<L a n L VHL + l)-A h (A h ±l). 

n L 

It is seen that the SM operator connects only compo- 
nents from one and the same basis level nL. Indeed, the 
operators L x and L y have no r-dependent part and so, 
due to orthogonality of the basis functions R n L(r), can- 
not connect the components with different nL. But the 
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latter can be done by the quadrupole operator r Y21. In 
this sense the SM operator is more selective than E21, 
though both operators generate transitions of the same 
multipolarity. K n = 1+ states involve both SM and E21 



modes and respond to both Mil and E21 external fields. 
The states are treated as magnetic dipole SM or electric 
quadrupole E21, depending on each of the two responses 
dominates. 
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FIGURE CAPTIONS 

Figure 1: Photoabsorption dipole cross section in 
light axially deformed Na clusters. Quadrupole and hex- 
adecapole deformation parameters are indicated in boxes. 
RPA results are given as vertical bars (in eVA 2 ) and 
as a strength function smoothed by the Lorentz weight 
with the averaging parameter 0.25 eV. Separate contri- 
butions to the strength function from the A/i =10 and 11 
dipole branches (the latter is twice stronger) are shown by 
dashed curves. The experimental data (triangles) from 
|'2.'"{l| are given for the comparison. 

Figure 2: Two-photon process: scheme of population 
of IRQM states A/i =20, 21 and 22 via the A/i = 10 (left) 
and 11 (right) branches of the dipole plasmon. 

Figure 3: Quadrupole strength distribution in light 
Na clusters. Lower panels: the unperturbed leh strength 
(without residual interaction). Upper panels: the RPA 
strength (with residual interaction). The results are 
given as bars (for every discrete RPA or leh state) and 
as smooth strength functions obtained by folding with 
a Lorentzian of width A =0.25 eV. The IRQM strength 
(enclosed by the circles at 0.5-1.5 eV) is very weak and 



so rescaled by the factor 10 2 . 

Figure 4: The electron level scheme for Na^", Na^ 
and Na^ 5 in the spherical limit (left) and at the equi- 
librium deformation (right). Occupied and unoccupied 
states are drawn by solid and dashed lines, respectively. 
The Fermi (HOMO) level is marked by index F. Arrows 
depict the possible low-energy hole-electron E2[i transi- 
tions. 

Figure 5: IRQM in light clusters. The plots exhibit 
quadrupole photoabsorption (uppermost panels, marked 
E2), scissors Ml photoabsorption (second line, marked 
Ml), and two-photon population of IRQM via A/i = 10 
(third line, marked TPP E10) and 11 (lowest panels, 
marked TPP Ell) dipole branches. IRQM are depicted 
by solid (A/i = 20), dashed (A/i = 21), and dotted 
(A/i = 22) curves. The strengths are smoothed by a 
Lorentzian with A = 0.1 eV. 

Figure 6: Quadrupole strengths for particular IRQM. 
The unperturbed leh (dashed curve) and RPA (solid 
curve) strengths are compared. The results are smoothed 
by a Lorentzian with the width A = 0.1 eV. 
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